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We extend the WKB method for the computation of cosmological perturbations during inflation 
beyond leading order and provide the power spectra of scalar and tensor perturbations to second 
order in the slow-roll parameters. Our method does not require that the slow-roll parameters 
be constant. Although leading and next-to- leading results in the slow-roll parameters depend on 
the approximation technique used in the computation, we flnd that the inflationary theoretical 
predictions obtained may reach the accuracy required by planned observations. In two technical 
appendices, we compare our techniques and results with previous flndings. 
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I. INTRODUCTION 

One of the main topics in theoretical cosmology is the 
computation of the spectra of perturbations during infla- 
tion (for an introduction and review see, e.g. Refs. 
which are then to be compared with the data for the Cos- 
mic Microwave Background Radiation (CMBR) obtained 
in present and future experiments It is in general 

impossible to obtain such spectra analytically, and there- 
fore approximate methods have been proposed, amongst 
which stands out the slow-roll approximation of Ref . 01 ■ 
This amounted to introducing a suitable hierarchy of pa- 
rameters, whose time-dependence was frozen to first or- 
der, and corresponding spectra and spectral indices were 
then determined uniquely. 

The task of extending standard slow-roll results be- 
yond first order is however highly non-trivial. As was 
noted in Ref. 0, one expects that the time-dependence 
of the parameters that describe the inflationary Universe 
could then appear in the order of the Bessel functions 
which describe the perturbations. Since this is, in gen- 
eral, inconsistent, other ways have been devised in order 
to by-pass this issue, such as the Green's function method 
of Refs. 0, Q ^nd the "uniform" approximation devel- 
oped in Refs. [111 [13. 

Our approach to this problem is to make use of an 
improvement of the WKB approximation (first applied 
to cosmological perturbations in Ref. introduced in 
Ref. 1^] . We followed the method of Ref. il3| in order to 
illustrate an improved (uniform) leading order approxi- 
mation and then introduced two expansions: the "adia- 
batic" expansion of Ref. and a perturbative Green's 
function expansion. In the present paper, we shall use 
our methods (in particular the adiabatic expansion of 
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Ref. ^3) ill conjunction with the slow-roll approxima- 
tion |j, [lj| in order to obtain higher order terms in the 
slow-roll parameters. Our method will not require that 
the slow-roll parameters be constant and extends the re- 
sults illustrated in Ref. Let us emphasize here the 
distinction between the slow-roll and WKB formalisms. 
The latter (also referred to as "semiclassical" or adiabatic 
expansion) is used to obtain approximate solutions to the 
general cosmological perturbation equations, whereas the 
former corresponds to an approximate solution for the 
homogeneous mode in which the inflaton kinetic term is 
neglected and one has a hierarchy of horizon flow func- 
tions (see Section Thus, they correspond to diverse, 
and independent, approximation schemes, even if, as is 
obvious, the validity (or not) of the slow-roll approxi- 
mation will influence the form of the approximate WKB 
solution. 

In Sectional we briefly recall the main definitions and 
standard slow- roll results. Our results to leading WKB 
and second slow-roll order, and next-to-leading WKB 
and first slow-roll order from Ref. are respectively 
reviewed in detail in Sections |llll and Hvl New results for 
the next-to-leading WKB and second slow-roll order are 
given in Section We finally comment on our work in 
Section IVII and compare with other results in the litera- 
ture in Appendices 1X1 and IbI 



II. BASICS 

Let us begin by recalling the Robertson- Walker metric 
in conformal time 77, 

ds2 = a^^ri) [~Aff + Ax ■ dx] , (1) 

where a is the scale factor of the Universe. Scalar (den- 
sity) and tensor (gravitational wave) fiuctuations are 
given on this background, respectively by /i = /is = ^ Q 
{Q is the Mukhanov variable _16J) and /i — /it = ah {h 
is the cimplitude of the two polarizations of gravitational 
waves [13, Il3), where the functions must satisfy the 
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one-dimensional Schrodinger-like equation 



together with the initial condition 
lim M(fc,ry)- 



(2) 



(3) 



In the above k is the wave-number, H — a'/a^ is the 
Hubble parameter, and 



where 



zs = a' 



' /H for scalar and z = zt 



(4) 



for 



tensor perturbations (0 is the homogenous inflaton and 
primes denote derivatives with respect to 77). The di- 
mensionless power spectra (PS henceforth) of scalar and 
tensor fluctuations are then given by 



fc3 

2 7r2 



ZT 



and the spectral indices and runnings 



dlnV, 



dln/c 



dlnP^ 
dlnfc 



" (dlnfc)2 ' = (dlnfc)2 

We also define the tensor-to-scalar ratio 



(5a) 



(5b) 



(5c) 



The evolution of the Universe is usually described by 
means of a set of fiow equations 0,^3 ■ The zero horizon 
fiow function is defined by 



H{N) 



where N is the number of e-folds, N = ln(a/ai) [where 
tti = a(r]i)] after the arbitrary initial time rji. The hier- 
archy of horizon flow functions (HFF henceforth) is then 
defined according to 



din I 



dN 



i > 



(7) 



and inflation takes place for ei < 1 (see Appendix fXlfor a 
comparison with different conventions) . We flnally recall 
that the frequencies can now be expressed in terms of 
ei, 62, and 63 only, as 



x2 7?2 



^=2-61 



(8) 



3 1 1 2 1 

a2i/2 =2-ei + 2^2-2^ie2 + -62 + -e2e3 



In general, fl'^{k, rj) vanishes for a certain time 77 = ?7o(A:) 
(the classical turning point). 

A. WKB Formalism 



Let us briefly summarize the essential formulae from 
our Ref. |^|. First of all, we apply the Langer transfor- 
mation 



a; = In 



Ha 



X^(l-ei)i/2e--/V (9) 



(5d) so that Eq. 10) becomes 



— +u''{x) 



(10) 



(6) where for scalar and tensor perturbations one now has 



ujI{x) 



LjJ^{x) = 



c,2x 



1 /3-eiV (3-2ei)e2 (l-2ei)e2e3 {l-4^^i)4 



(l-ei)2 4Vl-ei; 2(l-ei)2 2(l-ei)3 4(l-ei)4 



(11) 



,2 a; 



_ 1 (3_ei 
(l-ei)2 4 



1-ei 



ei £2 



6162 £3 , (2 + ei)eie 



2(l-ei)2 2(l~ei)3 4(l-ei)4 ' 
I 



and the corresponding turning points are mapped into the region to the right (left) of the turning point xo with 
X — xo{rj(),k). We then denote quantities evaluated in 
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the subscript I (II). For example, 



uji{x) = y/uj'^{x) , wii(x) = ^-uj'^{x) (12) 



and 



Ci(x) = / uji{y)dy , ^n(a;) = / ujn{y) dy . (13) 

J Xq J X 

The next-to-leading WKB solutions of Eq. (^UJ, in the 
adiabatic expansion and for one linear turning pointy is 
then given by a combination of Bessel functions and their 
derivatives as [l^lT^ 



C//±) {x)^[l + 01(1) {x)] {x) + 71(1) ix) u't^' {x) 



(14) 



U\^\x)=[l + c^n^i){x)\ u}f^(a;)+7ii(i)(x)u{i^^'(x) , 
where 



(15) 



The above expressions (|15|l satisfy an equation of motion 



dx-^ 



+ u/{x) — a{x) 



X{x) = , 



where the residual term is defined by 

3 {uj'f uf_ h_uP_ 

4 tj2 " 2^ ^ 36 ^ 



cr(a;) 



(16) 



(17) 



As we showed in Ref. ilj, such a quantity is expected 
to be small for the cases of interest, that is, in the suh- 
horizon limit 



aH 

in the super-horizon limit 
k _ 



-foo , 







(18) 



(19) 



and also at the turning points (zeros) x = a;o of the fre- 
quencies. In Eq. p7|) . we omitted the subscripts I and II 
for brevity and primes hereafter will denote derivatives 
with respect to x. 

The next-to-leading WKB PS in the adiabatic expan- 
sion are given by 



3 p2^ii,s 



neiml^\aHj {I - ei) uju^s 



16. 



(20a) 



i-Pi 



aH 



(1 - ei) cJiLT 



The spectral indices are also given by 



ns - 1 = 3 + 2 



d In /c d In A; 
AD 



(20b) 



TT-T = 3 + 2 



d^ii,T ^ '^■g(l)T 
d In fc d In fc 



and their runnings by 



as = 2 



ax = 2 



d 4ii,s , " 5(i)s 



{d\nkf {Alnkf 



(20c) 



d^ ryAD 
Q 5(i)T 



(dlnfc)^ (dlnfc)^ 



Finally, the tensor-to-scalar ratio takes the form 



?ii,s 



(20d) 



t^II.T 



where all quantities are evaluated in the super-horizon 
limit [x xq) and we used the results 



9{^s,t{^) = 2 I 011(1) (a;) - 7ii(i)(a;) 



UJll{x) 



- 01(1) (a;^) 



_tjjj(x)_" 
2tJii(x)_ 

(21a) 



S,T 



'^(i)(a;) 
7(1) (a;) 



2uj{x) J u}{y) 



7(i)(y) + c^(y) 7(1) (y) 

ct(2/) , 



Ay 



(21b) 



where it is understood that the integration must be per- 
formed from xq to X in region I and from x to xq in 
region II, and the indices S and T refer to the use of the 
corresponding frequencies. Let us note that if we take the 
above formulae with 5(\)g t(^) = 0, we obtain the lead- 
ing WKB results given in Ref. [13 for the PS, spectral 
indices and a-runnings. 



B. Slow- Roll Results 

Let us briefly recall the known results obtained from 
the slow-roll approximation 0, 0, Q| . To first order in 
the slow-roll parameters {7}, we have the PS 
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Trmpj 



l-2(C+l)ei-Ce2-(2ei + e2) In 



l-2(C+l)ei -2ei In 



(22a) 



where C = 7E+ln2-2w - 0.7296, 7^, being the Euler- 
Mascheroni constant, and aU quantities are evaluated at 
the Hubble crossing r;, , that is the moment of time when 
/c* = {aH)(N^) = (ai/)*. (The number is usually 
called "pivot scale" .) From Eqs. H22a() . we can also obtain 
the spectral indices and a-runnings, 

- 1 -2 ei - £2 , nT = -2 ei (22b) 



leading WKB approximation. 

We first set g^^{x) = in Eqs. plajl . (j20b|l . p1c|) 
and Ij20d|) . We then find it convenient to re-express all 
relevant quantities in terms of the conformal time 77. For 
this purpose, we employ the relation 



as = Q-T = 



(22c) 



on using respectively Eqs. (|5b|) and lf5c|) . From Eq. Ij5d|l . 
the tensor-to-scalar ratio becomes 



i? = 16ei 



1 



C £2 + £2 In ( — 



(22d) 



k rj 



k 



[1 + £1 £2 + O (e^] (23a) 



[f + £l£2+0(£^)] . (23b) 



III. LEADING WKB AND SECOND 
SLOW-ROLL ORDER 

We now consider the possibility of obtaining consis- We can now expand the frequencies to second order in 
tent second order results in the parameters £,;'s from the the HFF, 



UjUtj) = k^l]^ (l-2£i£2)- f^+3£i+4£? + ^£2+2£i£2 + j£^ + i£2£3 



'9 1 

w|(?7) = fc2^2 (l_2£i£2)- -+3£i-|-4£2--£i£2 



4 '2 



and the PS become 



H? , , ,3 2 / 8 1 19 n 19 1,1 \ 9, 



16 Hf . , ..c, 2 8 19 o 26 

r 



From now on, the subscript / will denote that the given form 
quantity is evaluated in the super-horizon limit. 



We must now compute the arguments of the exponen- > / ■ k) — /A^i ) — k^ — ^ ^-^ (26) 

tials in the above expressions, which are in general of the iiW0j?7fi J J^^ V vl) V ^ j 



(24) 



(25) 
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where the function A'^{r]) contains the HFF, but does not 
depend on k, and we have used 

dx=^ [l + e,{Tj)e2iv) + 0{e^)] . (27) 

T] 

Let us also note that, at the turning point rj = ijq, one 
has 



(28) 



It is now clear that, in order to obtain consistent results 
to second order in the slow-roll parameters, we must con- 
sider the time-dependence of the e^'s, and the function 
A'^(r]) may not be approximated by a constant. This 
does not allow us to perform the integral unless the scale 
factor a = a{r]) is given explicitly or, as we shall see, 
some further approximation is employed. 

In Ref. [l^, we proposed a procedure which will now 
be described in detail. Let us start with the general exact 
relation 



V2 



A{ri) ^ ( A{ri) - V^^(^) - v"^ 



+ 



In 



A{ri) + ^/A'^{ri) - rf 



A{t]) + ./A^(Jfj~n^ I 4:A{r]) 



dry (29) 



which holds for every function A^ {rj) . We can derive this 
relation by initially considering the integral in Eq. I|26l) 
with A constant (i.e. independent of rj, which is just 
power-law inflation). In this case, the integration can 
be performed and yields the first and second term in the 
r.h.s. of Eq. H29|l with the function A independent of r], 
whereas the third term vanishes identically. If we then 
reinstate the ry-dependence of A in the two non- vanishing 
terms and differentiate them with respect to 772, we ob- 
tain the original integrand in Eq. (|26|l plus another term. 



originating from A'^{ri). We then obtain the result H29|) 
by (formally) integrating the latter term and subtracting 
it from the previous two. 

One can also repeat the above procedure for the new 
integral in the r.h.s. of Eq. (|29|l . To do this, we define 



1 
V 



Biv) 



(30) 



and, after several simplifications, we can write the general 
relation as 



A{rj) + y/A^{r])-k^rf J 4 ^(77) 



8A{r,) 



lQriA^{q) SA{t])J 



with 



/ Ajrj) - y/A^{T]) - k^T]^ 

''U(77)-f 7^2(77) -fc2 ^2 J 8 7,^3(77) 



drj 

d77 ,(31) 



\^A{V)J U(7;) + v/I^M^IV 



-^Li2 



-Li2 



^(77) - ^A2(^) _ fc2 ^2 



A{r]) + ^^2(77) - fc2 772 
^A(^) 



(32) 



and Li2(z) is the dilogarithm function (see, e.g. Ref. |2g|). 



k=l ■^'^ 



In(l-z') , , , 

— ^ '-dz' . (33) 



On putting together all the pieces, we finally obtain 



rv2 



■ d rj 



\/A'^in) - fc2 772 — '- ^ JA'^M - fc2 



Yiv) 



Bin) 



AM 



8^(77) 



In 



A{r]) - y/A^Jr])~~Wrf 



A{t]) + y/A'^(r]) - fc2 772 



"n.) r^^+^)d77 

V2 

\A{7j) + ^A^r^) - k^ tj^ I SrjA^irj) 



A{r]) - ^A^r])~k^7f\ B^{t]) 



drj 



r 



(34) 



For the cases of interest to us [i.e. for the frequen- 
cies (|24|l ]. the fimctions A for scalar and tensor pertur- 



bations to second order in the HFF are given by 
3 . o 13 1 o 1 



42 



9 

- + 3.1 



. 2 13 1 2 

-e2 + 4ef + yeie2 + ^e2 



2^2 63 

(35a) 



- 3 ei -I- 4 ef + 4 ei 62 
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so that 

3 

i?s = 3ei £2 + 2 ^2 £3 , i?T = 3eie2- (35b) 
For the results to second order in the HFF, we can neglect 



the last two integrals in Eq. (PU, since B'{r]) = O (ef) 
and B^irf) = O i^f)- Thus our approximation consists in 
neglecting such terms and not in assuming the HFF be 
constant (or any specific functional form for them). We 
can finally write 6ls(?7o, m\ k) and <^ii,t(??o, Vf, k) as 



3 3 1n3 , ^ ln3 A o A 2 1 2 

6i.S - "2 + ^-+ln3eij + — e2,f + (li^3+-j elf + —elj 



11 In 3 In 3 



6 



2 
1 



24 ^ 3 



ei,f e2,f 



In 3 ln^3 
~6 4~ 



2 

48 



e2,f e3,f 



2 2 



ln3 
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/ln3 




j ei,f e2,f + 







In(-fcryf) 



fei.f e2,f + ^e2,f e3,f ) In^ (-km) 



(36a) 



3 3 ln3 
2 + ^ 
3 
2 



ei,f 



ln3 ci 



"i,f 



ln3+ - I e 



2 



4 1n3 ln^3 



In 3 



ln(-fc77f) - -ei,f e2,f In {-kru) 



(36b) 



where the arguments in have been omitted. On in- 
serting the above quantities in Eqs. H25|l . we can calculate 
the PS. 

In order to compare the PS from Eqs. H25I) . (|36al) and 
(|36bp . with the slow-roll expressions in Eqs. 122a|) . we 
need a relation between H{Nf) and H{Nf,), ei{Nf) and 
ei(^*)i and so on, to second order in the HFF (hereafter, 
quantities without a subscript will be evaluated at the 
Hubble crossing r/, corresponding to the pivot scale fc*). 
We expand the parameters ei^f to first order in AA^ = 
Nf — N^: in the numerators and to second order in the 
denominator of the scalar spectrum, 

~ 1 + e,+iA7V + l (e^^i -I- e,+i £,+2) AiV^ ,(37) 

and, analogously, we find 



iJ2 



1 - 2ei A7V- (eiea - 2ef) ATV^ . (38) 



We can eliminate rjf in the logarithms by expressing it in 
terms of l/(a H)[, 



ln(— fc 7/f ) 



In 



In 



k (1 -l-eij e2,f) (aff)* 



{aH){ (1 



+ eij A7V + ei,f 



In 



A7V + ei AiV + ei , (39) 

where in the first equahty we used Eq. H23a|l . in the sec- 
ond one the definition of the pivot scale, and, in the last 
relation, Eq. (|37|) . If we now use the above expressions we 
obtain the following PS, where the superscript stands 
for second slow-roll order and the subscript WKB for 
leading adiabatic order. 
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(^.WKB 



H 



1) £1 -i:'wKBe2 + 2D 



2 i^WKB 



9 



Dv 



12 



20 



ei £2 



[-2 ei - £2 + 2 (2 DwKB + 1) e? + (2 D^kb - 1) ei £2 + ^^wkb e 



2 WKB 



2 - n 

2 



" 24 " 18 
:e2e3] In 



£2 £3 



1 



(4£^ + 2£i£2 + £2-£2£3) In 



2 f _^ 



(40a) 



75 



(2) 

/l,WKB 



9 ^WKB 



TT m 



PI 



1 - 2 (Dv 



-^WKB 



2D^ 

^ '4£2-2£l£2) In' 



1) £1 



19 



2 ^WKB 2 Dv</KB g I £ 



£1 £2 + [-2 £1 + 2 (2 D„KB + 1) £? - 2 (D„KB + 1) £i £2] In (^^^ 



In the above, we have defined Dwkb = In 3 w —0.7653, 



which approximates the coefficient C in Eqs. H22a(l with 
an error of about 5%, and the factor Awkb = 18 /e^ ~ 
0.896 which gives an error of about 10% on the amph- 
tudes. One of the main results of this paper, aheady 
reported in Ref. is that the PS to second order in 



J 



r 



the slow-roll parameters do not depend on A7V, once the 
laborious dependence on the Hubble crossing has been 
worked out: this fact is in complete agreement with the 
constancy of the growing modes of C and h on large scales. 
The spectral indices, from Eqs. (|20b|l . and their runnings, 
from Eqs. H20c(l . are analogously given by 



(2) 

'S.WKB 



1 = -2 £1 - £2 - 2 £1 - (2 D„KB + 3) £1 £2 - -DwKB £2 £3 - 2 £1 £2 In 



(2) 

1 

T,WKB 



-2€i-2£2-2 (D„KB + 1) ei £2-2£i£2 In 



£2 £3 In 



(40b) 



J2) 



R 



(2) 



16 £1 



1 + Dv 



; £2 



-DwKB + g I £1 £2 



J2) 



"S.WKB - ~2£i£2 - £2£3 , "t.wkb 

From Eq. Ij20d|l the tensor-to-scalar ratio becomes 



= -2£l£2 



1 



I k \ 1 f k 

+ (£2 + £1 £2 + -D„KB £2 + -DwKB ^^2 £3) In ( — j + - (£2 + £2 £3) In' ( ^ 



(40c) 



1 2 _ ^ }_ 

2 ™ ~ 24 18 ' ^^^^ 



(40d) 



r 



We note that Eqs. iBUal) . I|40b|l and (plcl) agree with 
Ref. [m, to first order in the HFF. 

We can also compare the above results with others as 
we show in Appendix^ For example, one could approx- 
imate the HFF by their Taylor expansions around t/q ^-s 
in Ref. H, 0, ^ , and give the results in terms of the ei 's 
calculated at that point (see Table ITTll for results only on 



spectral indices). 
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IV. NEXT-TO-LEADING WKB AND FIRST 
SLOW-ROLL ORDER 

To first order in the HFF, the function A can be taken 
as constant. In this case, let us then introduce the more 
convenient notation 



£,±{xo,x;k) = ± / VT^^(y)di 

J X 



± 



1 ~2 



= ±A [atan± (e±) - e±] 
where we have used Ax = At/t^ and 

= ^+361 + ^62 



(41) 



(42a) 



3ei 



kijo = ~A, (42b) 
with ei and £2 constant as well. It is also useful to define 



atan_(-(a;) = arctanh(a;) 
atan_ [x) = arctan(a;) . 



(43a) 



(43b) 



where the plus (minus) sign corresponds to region II (I). 
With this notation, the expressions for 7(1) (77) and (p^i) (77) 
in Eqs. (|21b|) can be explicitly written as 



1 r 5T3ei 



A^ \ 24 72 e± [atan± (e±) - e±] 

(44) 

505 T 654 eg, + 153 9^ 5 [(102 t 90 6^,) atan± (9±) - (102 t 157 9|) Q±\ \ 
Mim - ^42 1^ ± 115201 10368 6^ [atan± (e±) - 



The limits of interest are then given by 



lim 7-(i)(»7)=0 



23 



lim 0_(i)(?7) - 

r)^-oo ^ ' 3150 



lim 7+(i)(r7) 

lim 0+(i)(?7) 
>)— »o 



1 



12 
181 
16800^2 



where we have used 



lim 9-(?7) =00 

rj — > — oc 



lim Q+{r]) = 1 , 



and 



lim tj„ + — 
7)^0- \ 2tJn 



A . 



(45a) 
(45b) 
(45c) 
(45d) 

(45e) 
(45f) 

(45g) 



We can now use the values in Eqs. (|42a|l to calculate the 
correction in the super- horizon limit (rj 0^) for 
scalar and tensor perturbations, obtaining 



,AD 



37 19 



.AD 



37 19 

324 ~ 243 



(46) 



Finally, from Eqs. (|20a|l and 1)46(1 . we can write the ex- 
pressions for the scalar and tensor PS, with the super- 
script for first slow-roll order and the subscript WKB* 
for next-to-leading adiabatic order, as 
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C,WKB* 



^2 



TT ei TO 



2 ^WKB* 
PI 



(1) 

/l.WKB* 



16 i/^ 



TT TO 



PI 



1 - 2 (DwKB. + 1) £1 - i^wKB. £2 - (2 £1 + £2) In 



1 - 2 (Dv 



(47a) 



1) £i-2ei ln( — 



where 



= — ln3 ~ —0.7302, which approxi- 
mates the coefficient C in Eqs. (|22a|l with an error of 
about 0.08%, and the factor Awkb. = 361/18 w 0.999 
which gives an error of about 0.1% on the amplitudes. 
We also obtain the standard slow-roll spectral indices and 
Qf-runnings, 



'S.WKB. 



1 = 



2 £1 £2 , "t,wkb. — 2 £1 (47b) ^Yie computation of the functions 



V. NEXT-TO-LEADING WKB AND SECOND 
SLOW-ROLL ORDER 

In order to give the results to next-to-leading WKB 
order and second slow-roll order, we should evaluate the 
corrections for scalar and tensor perturbations to 
second order in the e^'s. This, in turn, would require 



'±(1) 



and 



7±(i) 



a. 



(1) 

S.WKB* 



^T,WKB. 



, 



on respectively using Eqs. Ij20bp and (|20c(l . 
Eq. I|20dl) the tensor-to-scalar ratio becomes 



16 ei 



l + i^WK B. £2 



(47c) 
From 

(47d) 



Eqs. Ij21b|) . whose analytical calculation is extremely dif- 
ficult. We shall therefore try to give the expressions 
for the corrections by following a heuristic (and faster) 
method. 

Let us start from the first order slow-roll expressions 
given in Eqs. (|46l) and add all possible second order terms 
in the HFF. Some of such terms can then be (partially) 
fixed by making reasonable requirements, which will be 
explained below, so that the corrections to the PS read 



1 + fffi^s 



AD 
" 5(1)T 



361 4 / 1 

324 1 ^~57 



+ 



2 2 



£2,f £3,f 



4 262 
57^^^"^™-^ " 3249 



£l,f £2,f 



13 
1083 



t2.f 



4 2 

^ £i,f £2,t + ^ £2,f £3,f ) In (-fc m) 



(48) 



361 
324 



1 ei f 

57 ' 



361 



4 N 16 



£l,f £2,f - ^ £i,f £2,f In (-fc ?7f ) 



First of all, we have factorized the number 361/324 so as 
to recover the standard slow-roll amplitudes H22a(l with 
a very good accuracy. As already mentioned, the first 
order terms are the same as those in the results (|46|) of 
Section llVl evaluated however in the super-horizon limit. 
The form of the coefficients multiplying the second order 
monomials in the HFF have been partially determined 
by using the expressions given in Section IIVI [Eq. (|45a|) - 
( |45g| )] with Eqs. H35a() replacing Eqs. (|42a|) . By this pro- 
cedure, one expects to obtain results correct up to deriva- 
tives of , which can just contain mixed terms to second 
order. The numerical coefficients in front of t\ j and t\ j 
arc thus uniquely determined, whereas the coefficients 
multiplying the mixed terms ei.f e2.f and 62, t £3.f remain 



partially ambiguous. The last number in each square 
bracket arises directly from this procedure, whereas the 
other terms are chosen so as to match the "standard" 
dependence on Dwkb. in the PS, with 6s_t and ds left 
undetermined. Finally, we have also required that the PS 
do not depend on AA'', consistently with the constancy in 
time of the growing modes for h and C^. This requirement 
fixes the coefficients in front of the logarithms in Eqs. H48|) 
uniquely 23], and leads to spectral indices which do not 
depend on 63. t and ds. 

Proceeding then as in Section Hill from the above cor- 
rections we obtain the expressions for the scalar and ten- 
sor PS, 
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:,(2) 
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[-2 ei + 2 (2 1?wKB. + 1) £? - 2 (i^wKB. + 1) £1 £2] In 
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We also obtain the spectral indices Ij20bp and their runnings H20c|l . 



(49a) 
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SjWKB 
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(49b) 



a. 



(2) 

S,WKB* 



-2eie2 - £2 £3 , 



(2) 



-2£l£2 , 



(49c) 



and the tensor-to-scalar ratio becomes 



(2) 

WKB* £2 



2 WKB* 



253 
1083 



{bT - 6s) ^ 



71 



1083 



+ o -^WKB. 



£1 £2 



49 \ 



777 - ds — + —- £2 £3 
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57 722 / 



+ (£2 + £1 £2 + i^wKB. £2 + ^wKB. £2 £3) In (^Y^ + ^ (£2 + £2 £3) In^ (^Y^ I 



(49d) 



Let us end this Section with a few remarks. The un- 
determined coefficients 63, t and ds still appear in the PS 
(and their ratio R), but not in the spectral indices and 
runnings, which arc therefore uniquely specified by our 
(heuristic) procedure. In general, one may expect that 
the complete adiabatic corrections to second order con- 
tain HFF calculated both in the super-horizon limit and 
at the turning points (zeros) of the frequencies. How- 
ever, this is not relevant in the present context, since the 
possible numerical coefRcients that would multiply such 
terms can be eventually re-absorbed in the definitions of 
&S,T and ds- This could only be an issue (in particular for 



the spectral indices) if we considered third order terms 
(since it would then become important where the second 
order terms are evaluated). 



VI. CONCLUSIONS 

We have shown that the improved WKB treatment 
of cosmological perturbations agrees with the standard 
slow-roll approximation Q to within 0.1%, finally re- 
solving the issue of a 10 % error in the prediction of the 
amplitudes to lowest order which was raised in Ref. |llj . 
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The next issues are the inflationary predictions to sec- 
ond order in the slow-roU parameters. After the resuhs 
on the running of the power spectra [2l|, second order 
resuhs have been obtained using the Green's function 
method with the massless sohition in a de Sitter space- 
time |ii 0]. In a previous manuscript [l5l |. and in full 
detail here, we have employed the leading WKB approx- 
imation to obtain the scalar and tensor power spectra 
to second order in the slow-roll parameters. The key 
technical point is to use Eq. (|34l) . which allows one to 
express the power spectra on large scales as a function 
of the Hubble crossing quantities, but with no explicit 
dependence on AN. As one of our main results, we find 
that the polynomial structure of the power spectra in the 
Ei's obtained from the WKB method agrees with the one 
arising from the Green's function approach of Refs. [p,H. 

In Section EI as a further development of Ref. [13, we 
employ the next-to-leading WKB approximation to sec- 
ond order in the e^'s and, by requiring that the power 
spectra do not explicitly depend on AA^ (a property 
which instead was derived both in Sections IIIII and ITVjl . 
and using the expressions found in Section HVl we obtain 
unique expressions for the spectral indices and runnings. 
Our findings show that different ways of approximating 
cosmological perturbations show up at the leading or- 
der in the amplitude, but in the next-to-leading order in 
the derivatives of the power-spectra with respect to wave 
number k. The accuracy in the theoretical predictions on 
spectral indices to second order in the slow-roll param- 
eters, evaluated at fc = fc*, is now striking: the Green's 
function method coefficient C ~ —0.7296 is replaced by 
DwKB* — —0.7302 in the WKB method, leading to a 
precision of 1 part on 1000 in the predictions of the co- 
efficients of 0{ef) terms in the spectral indexes. 



APPENDIX A: SLOW-ROLL VS HORIZON FLOW 

In this Appendix, we compare the HFF we used in the 
text with some other (equivalent) hierarchies. 

In Table m we give the relations between the slow-roll 
parameters defined by some other authors and the HFF 
defined in Eq. 0. 



TABLE I: Relation between (some) slow-roll parameters and 
the HFF used here and in Ref. 11]. 



J. Martin et al. [14r 


E.D 


Stewart et al. \±,6}'° 


S. Habib et al. [8] 


S. 


Habib et al. [a, iO] 








S — €1 — ^ £2 




5l = ^ £2 — ei 


6 = I £2 £3 


d2 = -l 


ei £2 + 2 e? + i + i £2 ea 



In Ref. (]Jl ^ is used instead of C2. ''In Ref. |Jl ei, 5, and 



cf) 5/H 4> are used instead of e, 5i and 82 respectively. 

We also compare the Hubble-slow-roll (HSR) and the 
potential-slow-roll (PSR) parameters (see Ref. 22j) with 
the HFF. We recall that some authors redefine and 
as and and write only e, 77 and ^, which can 
be confusing. The relations between the ei's and HSR 
parameters are 

1 

en = ei , riH ^ ei~ -^£2 

(Al) 

Cff = ^1 - 2 '^l + 2 £2 £3 , 

and for the PSR parameters we have 



£2 



2(3-ei) 



6 £1 - f £2 - 2 -H I ei £2 - I ej - ^ £2 £3 



3-£i 



2£i--£2 



(2 ei - 6 - £2) [8 - 6 (4 -f 3 62) + £1 £2 (18 -h 6 £2 7 £3) - £2 £3 (3 -I- £2 + £3 + £4)] 



4(3- e,r 



(A2) 



2 1 

4£i - 3£i£2 + -£2£3 , 



where we have shown both the exact formulas and approximate relations to the order of interest. 

We complete this review with the first three exact connection formulae between the PSR and HSR parameters |2S 

\2 



£y _ (3 - tih) 



£ff (3 - eHY 



Vv 



3 (£g + m) -vh-^h 



(A3) 



9 {3eHVH +eH- ^hVh) - 3 {^VHen - vUl + <^h) + m<yH 



and so on with ajj the next HSR parameter. 



(3-£h)2 
I 



APPENDIX B: OTHER RESULTS IN THE 
LITERATURE 



able - in terms of the HFF in Tables and HTM 



We use Tableland Eqs. (|A1|I . to compare the results 
for spectral indices - and for PS and runnings when avail- 
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TABLE II: Results obtained with the Green's function method 7j suggested in Ref. ^Qj; C = ln2 + 7b - 2 ~ -0.7296 with 
7_B the Euler-Mascheroni constant. We display the original results for spectral indices, Q-runnings and tensor-to-scalar ratio 
evaluated at the pivot scale fc,. 

E.D. Stewart et al. \4], E.D. Stewart et al. [6], S.M. Leach et al. \7], A.R. Liddle et al. [22] . 
n = 7^ {l-2(C+l) .i-C.2+(2C^ + 2C+4-5) el 

+ [-2ei -62 + 2(2C + 1) el + (2C- 1) ei £2 + C ei - C ea ea] In (^) + \ (4 + 2 ei ea + ei - £2 £3) In" (^)} 
T'h = |l -2(C+1) ei + (2(72 +2C + 4 e? + (-C^ - 2 C -f- - 2) 6162 

-h[-2ei-h2(2C-hl) e? - 2 (C -Fl) ei ea] In (7^) + § (4e?-2eie2) ln'(r:)} 
ns - 1 = -2 ei - £2 - 2 e? - (2 C -)- 3) ei e2 - C £2 £3 , wt = -2 £1 - 2 £? - 2 (C + 1) £1 £2 
as = —2 £1 £2 — £2 £3 , aT — —2 £1 £2 

i?=16£i [l + C£2 + (C-4 + 5) £i£2+(iC'-ir + l) '^2^3] 



In Table im we summarize the results for scalar and ten- 
sor perturbations given in Ref. , as calculated with the 
method proposed in Ref. . We also refer to some other 
papers which give the same results, although expressed 
in different hierarchies. 

Let us then compare with the results obtained by 
the uniform approximation 8, 9, .10]. The integrand in 
Eq. (9) of Ref. can be obtained from our Eq. ^ by 
the change of variables 2: = In {—kr]), u = e~^/^ fj,, which 
of course differs from our Eq. jnj. Therefore, the inte- 
grands in our Eq. ifT^ and Eq. (9) of Ref. differ by 
0{ef) terms, once both are expressed in conformal time. 
We also note that in Refs. H, 0, 0, the authors just 
give expressions for the spectral indices which depend on 
k (i.e. without referring to a pivot scale k = k^). Fur- 
ther, it is not easy to see that they obtain the correct 
scaling in In (k/k^), which yields the a-runnings, and the 



Ei's are evaluated at the classical turning points for the 
frequencies. 



TABLE III: Results obtained with the uniform approxima- 
tion. 

S. Habib et al. \10] 

nS - 1 = -2£1 - £2 - 2£? - (f - 2tt) £1£2 " -Tt) £2 £3 
TIT = -2£l -2£f - - Itt) £1£2 



In Tables Uni we display the results for scalar and ten- 
sor perturbations given in Ref. |lOj |. As it was properly 
stated in Ref. jlflj , the results for scalar and tensor spec- 
tral indices given in Ref. H, 0] were not fully expanded 
to second order in the slow-roll parameters. 
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